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Introduction 

In  recent  years  several  methods  have  been  introduced  for  estimating  power  spectra 
with  considerably  greater  resolution  than  that  provided  by  the  conventional  periodogram  or 
the  Blackman-Tukey  windowed  Fourier  transform.  Included  among  such  techniques  are 
maximum  entropy  spectral  analysis  (MESA)  introduced  by  Burg  [1] , the  autoregressive 
model  (AR)  spectral  estimation  introduced  by  Parzen  [2] , and  the  method  of  maximum 
likelihood  as  demonstrated  by  Capon  [3] . Other  methods  offering  high  resolution,  which 
utilize  the  Fourier  transform,  are  described  by  Gerchburg  [4]  and  Papoulis  [5] . More 
recently,  another  spectral  estimator  has  been  introduced  by  Gray  [6] . 

While  none  of  these  spectral  estimators  have  been  thoroughly  investigated,  there  have 
been  a few  comparative  examinations  of  some  of  the  techniques  conducted  by  Lacoss  [ 7 ] , 
Ulrych  and  Bishop  [8] , and  Nuttall  [9] . Of  the  comparisons  investigated,  in  general,  superior 
results  are  achieved  by  using  the  MESA  method  and  the  Burg  technique  [10]  for  estimating 
filter  coefficients.  The  results  are  dramatic,  and  suggest  that  investigations  of  MESA  and 
other  high-resolution  techniques  be  continued.  Because  of  the  high  resolution  and  stability 
achieved  with  MESA  in  such  initial  investigations,  these  same  properties  are  investigated 
further  in  this  report  where  MESA  is  applied  to  the  analysis  of  simulated,  multichannel, 
spatial,  phased-array  data. 
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In  the  initial  paper  by  Burg  [1] , where  the  principle  of  MESA  is  first  suggested,  the  f 

prediction  error  filter  coefficients  (which  maximize  the  entropy)  are  specified  with  knowl- 
edge of  the  autocorrelation  coefficients.  However,  in  a second  MESA  paper,  Burg  [10] 

defined  the  prediction  error  filter  coefficients  as  a function  of  a set  of  uniformly  spaced  f 

data  samples  representative  of  the  function  of  interest.  In  addition.  Burg  simplified  the 
method  for  obtaining  the  filter  coefficients  with  use  of  Levinson’s  recursion  equations,  and 
also  noted  in  the  second  paper  that  the  mean  squared  prediction  errors  may  be  minimized 
in  both  the  forward  and  backward  directions.  These  suggestions  served  to  greatly  facilitate 
the  implementation  of  MESA  and  to  significantly  enhance  the  MESA  properties.  Taken 
together,  the  improvements  suggested  by  Burg  [10]  are  often  referred  to  as  the  “Burg 
technique.” 

The  concept  of  estimating  power  spectra  by  maximizing  entropy  appears  unique  in 
the  history  of  science,  yet  the  resultant  expression  for  power  spectra  is  identical  to  the 
representation  of  the  all-pole  method  or  autoregressive  model  (AR)  introduced  by  Parzen 
[2] . In  fact,  van  den  Bos  [11]  and  Kaveh  and  Cooper  [12]  have  noted  that  MESA, as  out- 
lined by  Burg  [1] , is  equivalent  to  the  AR  method  as  described  by  Parzen.  Therefore,  it  is 
of  consequence  to  note  that  the  different  spectral  estimates  sometimes  predicted  with  the 
two  spectral  estimation  methods  are  not  due  to  inherent  model  differences,  but  rather  to 
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the  different  methods  used  for  evaluating  the  corresponding  filter  coefficients.  With  this 
realization,  Ulrych  end  Bishop  (81  conducted  a comparative  analysis  of  the  Burg  and  Yule- 
Walker  113,141  techniques  for  evaluating  MESA  and  AR  filter  coefficients.  In  a comparison 
of  spectral  estimation  of  harmonic  components,  Ulrych  and  Bishop  showed  that  the  MESA- 
Burg  technique  provided  significantly  greater  resolution  than  did  the  AR-Yule-Walker  tech- 
nique. They  noted  that  the  resolution  differences  are  not  surprising  when  it  is  realized  that 
the  Yule-Walker  technique  has  assumed  a zero  extension  of  the  data  samples,  whereas  the 
Burg  technique  contains  no  assumptions  concerning  the  nonsampled  data  field. 

Since  the  MESA  and  AR  methods  are  most  significant  when  processing  short  daU  seU, 
it  is  natural  to  use  such  methods  for  processing  data  collected  with  multielement  electro- 
magnetic antennas  or  acoustical  arrays.  Such  methods  may  well  make  it  possible  to  achieve 
high  resolution  by  using  unusually  short  antenna  arrays.  Consequently,  in  this  report  the 
resolution  property  of  the  MESA-Burg  technique  is  examined  as  a function  of  input-daU 
signal-to-noise  ratios  (S/t?),  number  of  antenna  elements,  numbers  of  signals,  incident  signal 
angle,  relative  signal  amplitudes,  and  relative  signal  phase. 


THEORY 
Power  Spectra 

The  maximum  entropy  power  spectra  (MESA),  which  was  introduced  by  Burg  (1) , 
has  been  derived  by  Parzen  [21  using  statistical  methods  and  by  van  den  Bos  (11)  using  an 
all-pole  model  representation.  To  further  the  understanding  and  acceptance  of  MESA, 
another  derivation  suggested  by  Blizard  (15)  is  presented  based  upon  discrete  convolution 
filtering  and  minimization  of  the  mean  squared  error. 

Consider  the  one-step  discrete  convolution  prediction  filter  described  by  Levinson 
(16)  as  follows: 


(1) 


where  x is  the  prediction  at  time  t of  the  function  which  is  sampled  at  time  intervals 
n(At).  T^he  N prediction  coefficients  are  given  by  o^,  a^,  ■ ■ ■ The  error  of  the  one-step 
prediction  is  e^,  given  by 


The  filtering  and  error  analysis  represented  by  Eqs.  (1)  and  (2)  is  illustrated  in  Fig.  1,  where 
filter  coefficients  are  multiplied  by  values  of  x^  sampled  at  time  intervals  n(At),  and  the 
resultant  multiplications  are  summed  to  obtain  the  prediction  The  predicted  signal 
x^  and  the  actual  signal  x^  are  subsequently  subtracted  to  obtain  the  prediction  error  e^. 


■'t  • 
N 

- S 


(2) 
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Fig.  1 — Prediction  error  filter 


For  convenience,  the  prediction  error  filter  7 is  introduced  as  follows: 

N 

^ n • 

n“0 

where,  in  comparing  Eq.  (2)  and  (3),  it  is  observed  that 

\ ^ ^ -“fi*  ” ^ 1- 

Equation  (3)  is  transformed  to  frequency  space  with  the  Fourier  transform  to  obtain  the 
following  equation: 

s 

= Xicu)  X)  (4) 

fi”0 

where  the  Fourier  transform  of  a function  delayed  n(At)  units  is  the  exponential 

gi<jn(Af) 

multiplied  by  the  transformed  function.  The  power  spectra  density  function  Pfcj)  is  defined 
as 

/>(cj)  = X^CwVlF, 

where  W is  the  bandwidth  defined  by  the  sampling  interval  At;  i.e.. 


If  the  signal  distribution  function  X(co)  is  given  by  Eq.  (4),  the  power  spectra  density  func- 
tion becomes 
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P(u))  =• 


f^(co)/W' 


N 

I 

n-O 


in{^t)  2 


P(u;)  = 


p^/w 


1 + 


t y„' 


(5) 


(6) 


,iij  n(«^0  1 2 


n*l 


where  the  error  power  {^(u;)  is  represented  by  A requirement  that  the  spectral  error 
power  be  a minimum  results  in  Pjy  being  independent  of  frequency  as  follows: 


doj 


= 0. 


^^(cj)  = constant, 


Pjy  = constant. 


(7) 


Then  if  is  a constant,  the  prediction  error  filter  7 is  a whitening  filter,  and  is  also  the 
mean  of  the  total  squared  error  as  follows: 

Ps  = 

Pn  (8) 


The  power  spectra  density  function  P(cj)  given  by  Eq.  (6)  is  the  same  MESA  equation  intro- 
duced by  Burg  ( 1 ) and  later  derived  with  detailed  steps  by  Barnard  [17] . 

The  power  spectra  density  (Eq.  (6))  may  also  be  expressed  as  a wavenumber  power 
spectra  density  by 


P(k) 


N 


n-1 


(9) 


where  the  time  (f)  and  frequency  (/)  variables  have  been  transformed  to  space  and  wave- 
number  variables  by  using  the  following  relations: 


w = fee.  At  = Ax/c, 
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k “ (2»r/X)  cos  0, 

X » signal  wavelength, 

0 = signal  angle  of  incidence, 
c * signal  velocity, 

K - 2it/\. 


Prediction  Error  Filter 

Utilization  of  the  MESA  power  spectra  equation  (Eqs.  (6)  and  (9))  requires  that  the 
pretliction  error  coefficients  and  the  mean  error  power  be  known.  The.se  unknown 
parameters  may  be  specified  by  minimizing  the  average  time-dependent  pretliction  error 
power  e2.  The  resulting  N 1 equations,  which  are  derived  in  Appendix  A,  are  presented  in 
a matrix  formulation  as 


'"o 

'•2 

'■3  • • 

• ' 

>0 

V 

• 

'•l 

'•0 

rg  . . 

• '“n  1 

0 

• 

'•2 

• 

'•l 

'•0 

Tj  . . 

• ''at  2 

• 

>2 

= 

0 

* 

'•/V 

v-1 

'‘n  2 

'‘n-3  • 

• '‘0 

'>'n 

0 

(10) 


where  it  is  known  that  7o  = 1 and  it  is  assumed  that  the  autocorrelation  coefficients  r. 
(with  lag  (At)  are  known  for  N lags. 

The  autocorrelation  coefficients  have  the  following  definition: 


r 


n 


* Um  1 y y* 

M-»oo  2M  + 1 k n ‘ 


(11) 


But  for  a finite  data  set,  the  autocorrelation  coefficients  may  be  computed  by  approximat- 
ing Eq.  (11)  with  a finite  summation  over  M data  samples. 
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For  large-sized  filters  {N  large),  solution  of  the  N 1 equations  given  by  Eq.  (10) 
becomes  very  tedious.  Fortunately,  Burg  (10)  demonstrated  a more  expeditious  method 
for  specifying  the  unknown  prediction  error  coefficients  which  appear  in  Eqs.  (6)  and  (9). 
He  noted  that  the  unknown  parameters  and  may  be  evaluated  with  repeated  use  of 
Levinson’s  recursion  relations 


for 


(12) 


(13) 


N > 1, 


and  with  knowledge  of  which  is  shown  in  Appendix  B to  have  the  following 

representation: 


'N*l 


(14) 


where  the  forward  prediction  error  is  and  the  backward  prediction  error  is  The  three 
equations,  Eq.  (12),  (13),  and  (14),  comprise  the  Burg  technique  as  originally  demonstrated 
by  Burg  (10)  and  later  generalized  in  detail  by  Anderson  (18)  and  Barnard  (17) . 


The  remainder  of  this  report  is  concerned  with  the  properties  exhibited  by  the  MESA 
wavenumber  power  spectra  equation  (Eq.  (9))  when  evaluated  using  the  Burg  technique 
given  by  Eqs.  (12-14). 


Signal  Simulation 


Resolution  properties  of  the  MESA-Burg  technique  are  examined  by  using  simulated 
antenna  data.  Input  signals  to  a linear  (line),  multielement,  phased-array  antenna  are  assumed 
to  be  plane  waves  mixed  with  white  Gaussian  noise.  The  multichannel  signals  are  pre- 
processed  with  narrowband  filtering  and  channel  delays  which  serve  to  “direct”  (or  steer) 
the  antenna  in  the  direction  of  the  incident  plane-wave  signal.  The  preprocessing  methods 
are  illustrated  in  Fig.  2 where  the  nth  channel  is  depicted  as  delayed  (n  — 1)  A(  seconds  for 

At  = Ax  sin  (0)/c, 
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2 — Preproceuing  of  ■ linear  phased  array 


where  c is  the  velocity  of  the  incident  signal.  Ax  is  the  antenna  element  spacing,  and  0 is 
the  look  angle  of  the  steered  array  measured  relative  to  the  normal  to  the  array. 


The  signal  Jr^(0)  incident  to  the  nth  antenna  element  is  represented  as 

x^(0)  = exp(/J2„(D)J  + Q^exp{i2TrqJ,  (15) 


where  A is  the  signal  amplitude,  is  the  signal  total  phase,  and  is  a random  number 

representative  of  Gaussian  noise.  The  amplitudes  A and  Q are  relative,  and  are  determined 
by  a Gaussian  distribution  with  variance  and  a specified  signal-to-noise  ratio  as  follows: 


(16) 

(17) 


where  S/t)  is  the  input  signal-to-noise  power  ratio,  is  the  average  noise  power,  and  q is 
a random  number  between  0 and  1.0  obtained  by  using  a “white”  random  number  " 
generator. 


The  signal  phase  i2^(0)  has  three  components  as  follows: 

njO)  = 27ri[n -l](Ax/Xl  (sin  (0)- sin  (0^)1 -0,  (18) 

where  0 is  the  array  look  angle,  0^  is  the  angle  of  the  incident  signal  relative  to  the  normal 
to  the  array,  0 is  the  incident  signal  initial  phase  with  values  between  0 and  2n,  and  X is 
the  signal  wavelength  given  by 

X = c/f. 

For  all  signals  analyzed  in  this  report,  the  ratio  AxA  has  the  value  0.5;  i.e., 

AxA  = 0.5. 
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In  fonventionaJ  beamforminK  all  jV  a’.'*-**nna  elements  are  summed  such  that  the  total  \(0 ) is 

S 

X(tf)  « ^ expliS2^((/)l  ♦ exp  (i27r</^)[,  (19) 

n*l 

and  the  conventional  antenna  power  pattern  is  computed  in  decibels  as  follows: 

dB  = 10  lot?  (X^  («)/X^  (OJ).  (20) 

However,  the  MESA  technique  requires  multichannel  data  which  are  given  by  Kqs.  (15) 
and  (18). 

Conventional  antenna  patterns  are  compari'd  with  MESA  patterns  whenever  such  com- 
parisons are  considered  worthwhile.  It  is  specifically  noted  that  the  S/r}  is  defined  at  the 
antenna  element  and  is  the  same  for  each  antenna  element.  Consequently,  the  S/tj  does  not 
include  the  conventional  antenna  gain  factor. 


'4 


Power  Spectra  Peaks 


The  wavenumber  power  spectrum  given  by  Eq.  (9)  may  be  expressed  by 


. (21) 


i r 


whore 


i'(fe) 


N 

I + ^ ^ ^ ^(fcn(Ax) 
n«l 


7o  = 10. 
k = K sin  0, 
K = 2n/X. 


m) 


_ ^ifen(Ax) 

* fl 


1 

1 


Spectral  peaks  of  P(k),  which  occur  at  minimal  values  of  P(fe),  may  be  located  by  solving  for 
the  roots  of  the  function  ['(r)  in  the  complex  z-plane  as  follows. 
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(’onsid»*r  a polynomial  I'li)  dt*fimil  as 

V 

whort* 

2 = I2  , 

0 = h n ( Ajf),  (23) 

C>  = itn  sin  0 for  (Ax)  » X/2. 

The  polynomial  I’lz)  (of  order  .V)  may  also  be  expressed  as  a function  of  the  N complex 
roots  as  follows: 


1(2 ) » l(; -2j ) (2 • (2 -2^)1  , (24) 

where  z ^ , Z2,  ■ ■ ■ ■Zy  are  the  .V  roots  locat«*d  in  the  complex  z-plane.  The  significance  of  I'lz) 
becomes  obvious  from  the  observation  that  the  function  of  interest  l'(fe)  is  actually  the 
|3olynomial  I'fz)  on  the  unit  circle,  i.e., 

r(fe)  = 1(2)  for  Iz  I = 1,0. 

Consequently,  r(fe)  has  a polynomial  repn'sentation  in  the  z-plane  as  follows: 

[’(fc)  = l(e"^  -z,)(e"^  - Z2)...{e‘'^  (25) 

where  it  is  recalled  that  0 and  k are  functionally  r*>lated  according  to  Eq.  (23).  The  product 
terms  of  l'(fe)  in  P^q.  (25)  are  the  complex  vectors  {e“^  which  are  shown  constructed  in 
the  z-plane  diagram  of  Fig,  3.  In  accordance  with  Eq.  (25),  the  minimal  values  of  l'(fe) 
occur  at  minimal  values  of  the  complex  vectors  (e"^  — z ).  And  as  observed  in  Fig.  3 the 
complex  vectors  (e"^  — z^)  are  minimal  for  0 corresponc/ing  to  6 , the  argument  of  the  com- 
plex roots  zy,  i.e., 

<fi  = (pj  where  (pj  = arg  (Zj). 

Consequently,  stronger  spectral  peaks,  which  occur  for  smaller  values  of  l'(fc),  occur 
for  roots  zj  lying  closer  to  the  unit  circle.  .And  as  noted  by  Burg  (in  private  communication), 
all  roots  of  Eq.  (22)  lie  outside  the  unit  circle. 
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APPLICATIONS 
MESA  Examples 

Several  examples  of  MESA  antenna  patterns  are  computed  in  order  to  demonstrate 
MESA  characteristics.  Antenna  patterns  are  computed  for  linear  antennas  containing  8 or  16 
elements  spaced  at  half-waveicngth  intervals,  and  for  various  signal-to-noise  ratios.  Snapshot- 
type  antenna  patterns  are  presented  to  demonstrate  the  explicit  MESA  characteristics.  A 
MESA  snapshot  pattern  is  a MESA  antenna  pattern  computed  by  using  N discrete  data  points 
representing  spatial  data  collected  at  one  instant  of  time.  In  some  examples  several  MESA 
snapshot  patterns  may  Vie  averaged  to  obtain  one  single  representative  stable  antenna  pattern. 
Some  type  of  averaging  is  usually  recommended  when  MESA  is  applieKl  to  actual  data. 
Usually,  it  is  better  to  use  a time  average  of  either  the  input  autocorrelation  matrix,  the  pre- 
diction error,  or  the  prediction  filter  coefficients.  However,  in  this  demonstration  of  MESA 
characteristics  it  is  convenient  to  obtain  pattern  stability  with  the  averaging  of  snapshot 
patterns  (or  postprocessor  integration). 


The  MESA  algorithm,  used  to  compute  the  examples  in  this  report,  utilizes  complex 
number  arithmetic  throughout,  so  that  the  input-signal  and  noise  data  may  have  a complex 
representation.  As  a result,  the  prediction  filter  coefficients  computed  by  the  algorithm  are 
complex  coefficients. 


Each  example  MESA  antenna  pattern  is  computed  for  a noise  field  identified  by  a 
seed  number  required  of  the  random  number  generator.  The  exhibited  antenna  patterns 
may  be  duplicated  only  with  use  of  the  particular  seed  number  (IR)  that  is  identified  in 
each  example. 

Single  Signal 


A MESA  snapshot  antenna  pattern  is  shown  in  Fig.  4 for  a single  signal  incident  to  an 
8-element  antenna  at  0°  (broadside  to  the  antenna)  with  a S/r)  of  0 dB.  In  this  particular 
snapshot  pattern,  which  is  computed  by  using  a 6-point  prediction  filter  (N=6),  the  signal 
peak  is  about  10.5  dB  above  the  largest  noise  peak.  There  are  N peaks  in  a MESA  snapshot 
which  is  computed  for  an  N-point  prediction  filter,  and  in  Fig.  4 there  are  five  prominent 
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(H-aks  aiui  ono  small  poak  loi'att*Li  at  90“ . 'I’he  noiso  peaks  would  be  distributed  differently 
in  a MKSA  snapshot  eomputed  for  a subsecpient  noise  field.  'I'he  noise  field  ustnl  in  comput 
init  the  MKSA  snapshot  pattern  consists  of  einht  random  (white)  complex  numbers  repre- 
SL'ntative  of  receiver  noise. 

The  same  noise  data  of  Fin.  -I  summed  with  a signal  at  O"  having  a S/i)  of  10  dB  to 
obtain  the  MKSA  snapshot  pattern  shown  in  Fig.  5.  Since  the  noise  data  are  identical  in 
Figs.  4 ami  f),  the  noisi'  peaks  have  the  same  location  in  both  examples,  hut  the  noise  peaks 
of  Fig.  5 are  further  suppn'sst'd  (both  snapshot  patterns  are  normali^tHl  at  the  signal  peak). 

A conventional  antenna  pattern  is  also  shown  in  Fig.  5 for  comparison  purposes,  where 
the  conventional  pattern  is  ohservwl  to  have  a heamwidth  of  about  15“  at  the  half-power 
(iUIH)  points  on  the  main  lobe.  The  width  of  the  MKSA  signal  peak  in  Fig.  5 is  so  narrow  as 
not  to  he  measurable  on  the  plotteil  pattern.  Since  points  are  computinl  at  1“  intervals, 
it  can  only  he  stated  that  the  MKSA  signal  peak  is  down  by  10  dB  at  the  l“  intervals  on 
either  side  of  the  peak  oi’curring  at  O".  In  the  particular  MKSA  snapshot  illustrated  in  Fig.  5 
there  is  also  considerable  reduction  in  side-peak  levels  in  comparison  with  the  conventional 
pattern.  However,  side-peak  levels  vary  significantly  in  subs»Hiuent  MKSA  snapshots,  (’onsi.s- 
tenl  side-peak  levels  are  possible  only  with  some  type  of  preprocessor  or  postprocessor 
averaging. 

'/’ICO  Cltm'ly  Spaa'd  Signals 

Two  signals  incident  at  angles  of  4“  and  4“  are  readily  identified  in  the  MKSA  snap- 
shot shown  in  Fig.  6 for  a complex  noise  and  a S/rj  of  13  dB,  each  signal,  each  element.  The 
MKSA  snapshot  is  computed  for  an  8-element  antenna  and  a 6-point  prt'diction  filter,  which 
isdenotiHl  by  the  six  peaks  a)>|>earing  in  Fig.  6.  Since  particular  snapshot  patterns  are  not 
II  .'essarily  repetitive  over  a finite  time  interval,  a more  repivsentative  MKSA  antenna  fiat- 
lern  is  obtained  by  averaging  20  such  snapshot  patterns,  as  shown  in  Fig.  7. 
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Fig.  6 — Two  signals  at  —4*  and  4° 
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FiK.  7 — Two  signals  at  —4°  and  4°  (summed  patterns) 


While  the  two  signals  are  resolvt>d  in  the  averaged  MKSA  pattern  of  Fig.  7,  the  resolu- 
tion of  the  averaged  pattern  is  not  as  favorable  as  the  particular  MKSA  sna[Jshot  illustratt*d 
in  Fig.  6.  However,  the  two  signals  are  not  resolved  by  the  conventional  antenna  pattern 
(also  shown  in  Fig.  7),  which  is  computed  as  the  sum  of  two  conventional  patterns  obtaini*d 
by  “steering”  the  array  in  the  directions  of  the  two  signals.  It  is  known  that  the  resolution 
of  a conventional  antenna  pattern  is  given  apjiroximately  by  tbe  function  R as  follows: 

R ^ arc  sin  (A/L), 

where  X is  the  signal  wavelength  and  L is  the  array  length.  For  an  8-element  antenna  having 
half-wavelength  element  spacing, 


R ^ are  sin  (2/7) 

16.6  degrees. 

For  the  parameters  illustrated  in  Fig.  7,  the  resolution  of  MESA  is  at  least  a factor  of  two 
letter  than  the  conventional  pattern. 

Thrt'e  Closely  Spaced  Signals 

A MESA  snapshot  is  shown  in  Fig.  8,  computed  for  three  signals  incident  at  angles  of 
- 7°,  0°,  and  7°,  with  each  signal  having  a S/p  of  40  dB,  each  element.  The  MESA  snapshot, 
computed  for  a 6-point  filter  and  omnidirectional  noise,  has  four  equally  strong  peaks  and 
two  low-level  side  peaks.  It  is  apparent  that  the  central  signal  at  0°  is  repre.sented  by  two 
split  peaks  on  either  side  of  the  true  signal  location.  The  MESA  snapshot  pattern  depicts 
"peak  splitting,”  which  is  sometimes  observed  in  MESA  snapshots  at  the  larger  filter  sixes, 
but  is  encountered  less  frequently  in  averaged  MESA  antenna  patterns. 
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Pig.  8 — Splitting  with  three  signals 


For  example,  an  average  of  20  MESA  snapshots,  including  the  one  shown  in  Fig.  8 
(since  the  same  seed  number  is  used),  is  shown  in  Fig.  9 to  illustrate  that  pattern  averaging 
serves  to  enhance  MESA  stability.  In  comparing  Figs.  8 and  9,  it  is  observed  that  even 
though  stability  is  enhanced  by  averaging,  peak  widths  are  increastni. 

Further  examples  of  three  resolved  signals  are  shown  in  Figs.  10  and  11,  where  the 
effect  of  increased  S/rj  is  depicted.  In  Fig.  10  the  noise  field  is  omnidirectional  and  the  three 
peaks  are  predicted  within  0.5°  of  the  correct  incident  angles.  Also,  MESA  patterns  com- 
puted with  omnidirectional  noise  have  an  observed  symmetry  in  peak  location,  whereas 
MESA  patterns  computed  for  complex  noise  distributions  do  not  exhibit  such  symmetry  as 
demonstrated  by  the  snapshot  patterns  in  Fig.  11. 

The  MESA  snapshot  patterns  in  Figs.  10  and  11  do  indicate  that  resolution  is  improved 
and  side-peak  levels  are  reduced  with  increased  S/tj.  These  benefits  are  demonstrated  in  Figs. 
10  and  11  for  a S/rj  increase  of  20  dB.  The  two  examples  also  illustrate  that  omnidirectional 
noise  does  not  affect  the  accuracy  of  peak  location,  but  complex  noise  may  cause  such  in- 
accuracies or  distortions.  It  appears  that  complex  noise  (receiver  noise)  is  a more  troublesome 
noise  than  the  spatial  omnidirectional  noise  field. 

Pattern  Averaging 

While  it  has  been  demonstrated  that  pattern  averaging  provides  stability  to  Ml^SA 
antenna  patterns,  averaging  also  serves  to  broaden  peak  widths  in  comparison  to  MESA 
snapshot  patterns.  The  averaged  pattern  of  Fig.  9 demonstrated  MESA  stability  for  omni- 
directional noise. 

Another  benefit  of  pattern  averaging  for  omnidirectional  noise  is  illustrated  in  Fig.  12, 
where  side-peak  levels  have  been  substantially  reduced.  In  the  example,  a MESA  snapshot 
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pattern  is  shown  for  one  signal  incident  at  —3°  to  an  8-element  antenna.  From  the  number 
of  peaks  in  the  snapshot  it  is  noted  that  a 5 point  prediction  filter  is  used.  An  average  of 
20  such  snapshot  patterns  indicates  a substantial  reduction  in  side-peak  levels  that  is  much 
greater  than  anticipated  from  the  gain  due  only  to  a uniform  distribution  of  side  peaks  in  the 
20  MESA  snapshots.  For  in  an  ordinary  summation  of  snapshots  such  a gain  is  only 

10  log  (A/), 

where  M is  the  number  of  patterns  summed.  Consequently,  it  is  anticipated  that  side-peak 
levels  may  be  reduced  by  only  13  dB  in  an  average  of  20  MESA  snapshots.  No  explanation  is 
suggested  for  the  gain  that  is  apparent  in  F'ig.  12,  but  if  the  apparent  gain  is  real,  then  signal 
detection  may  be  improved  substantially  with  pattern  averaging. 

The  effect  of  averaging  MESA  snapshots  computed  for  complex  noise  is  demonstrated 
in  Figs.  13  and  14.  A MESA  snapshot  of  one  signal  incident  at  —3°  is  illustrated  for  a 7-point 
filter  in  Fig.  13  where  the  S/tj  is  0 dB.  It  is  apparent  from  the  averaged  MESA  pattern  shown 
in  Fig.  14  that  signal  detection  in  complex  noisc'  is  improved  with  averaged  MESA  patterns, 
but  the  amount  of  improvement  is  not  demonstrated  by  a single  MESA  pattern.  However, 
detection  performance  may  l)e  measured  with  a large  collection  of  such  patterns  by  con- 
structing a set  of  ROC  (receiver  operating  characteristic)  curves  that  illustrate  the  probability 
of  detection  for  specific  false  alarm  rates.  In  the  particular  example  of  Fig.  14,  there’ is  an 
apparent  gain  of  about  5 dB  in  S/rj  level  for  a zero  false  alarm  raU>.  The  gain  appears  even 
more  substantial  in  some  following  examples. 

The  S/t}  gain  that  is  apparent  with  the  previous  examples  of  averaged  MESA  patterns 
appears  to  indicate  that  the  gain  is  substantially  greater  with  MESA  patterns  computed  for 
omnidirectional  noise.  One  obvious  reason  for  this  discrepancy  is  observed  in  Figs.  10  and 
11,  where  complex  noise  affects  the  accuracy  of  the  location  of  signal  peaks,  thereby  reduc- 
ing the  effect  of  pattern  summation.  There  is  little  or  no  observable  effect  upon  peak  loca- 
tion with  omnidirectional  noise.  Consequently,  the  signal  peaks  summed  in  patterns  for 
omnidirectional  noise  remain  sharply  defined,  whereas  the  signal  peaks  summed  in  patterns 
for  complex  noise  are  distorted  and  broadened  and  show  only  a moderate  S/rj  gain  due  to 
pattern  averaging. 

To  recount  the  benefits  attributable  to  pattern  averaging,  the  examples  presenU'd 
indicate  that  averaging  of  a sufficient  number  (maybe  20  or  so)  of  MESA  snapshot  antenna 
patterns  provides  pattern  stability,  helps  to  eliminate  peak  splitting,  and  may  provide  a 
significant  gain  in  signal  detection. 

Mixed  Signals 

A MESA  snapshot  antenna  pattern  is  shown  in  Fig.  15  for  a collection  of  five  signals  in 
complex  noise  having  S/r/’s  of  0, 10,  and  20  dB  as  indicated  in  the  illustration.  The  example  of 
Fig.  15  is  computed  for  an  antenna  with  16  elements  and  a 13-point  prediction  filter.  Only 
two  of  the  signals  are  clearly  identified  in  the  MESA  snapshot  of  Fig.  15,  but  all  five  signals 
are  readily  and  accurately  detected  in  an  average  of  20  such  MESA  snapshots  as  shown  in 
Fig.  16.  Signal  detection  is  considerably  improved  in  the  averaged  MESA  pattern  which  con- 
sists of  signals  mixed  with  complex  noise.  It  is  apparent  that  the  particular  set  of  complex 
noise  samples  (IR  = 5)  has  not  adversely  affected  signal-peak  location  accuracy  as  was  noted 
in  the  previous  examples  of  Figs.  11  and  14. 
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Fig.  13  — Signal  at  — 3°  (snapshot) 
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Fig.  15  — Five  signals  (snapshot) 
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Relative  Signal  Phase 

It  has  htHMi  nottti  by  FouKen*,  et  al.  1 19)  that  signals  having  initial  phases  other  than 
zero  are  not  accurately  predictwl  hy  MKSA.  In  particular,  Fougere  noted  that  signals  with 
an  initial  phase  of  7r/-l  are  most  troublesome  for  MKSA,  as  peak  splitting  is  apparently  most 
s**vere  for  this  particular  phase. 

In  Fig.  17,  an  average  of  20  MKSA  snapshot  patterns  is  shown  for  the  same  five  signals  ' 

and  noise  and  the  same  antenna  and  filter  size  depicted  in  Fig.  16.  However,  the  signal  in- 
cident at  60°  has  a phase  of  rr/  l ( 45°)  relative  to  the  other  four  signals,  all  of  which  have 

zero  initial  phase.  In  comparing  Figs.  16  and  17  it  is  noted  that  the  two  signals  locaU*d  at  ' 

30°  and  60  are  considerably  broadeniKl  hy  the  [ihase  shift,  which  suggests  that  some  of  the  ^ 

summed  MKSA  snajishots  either  contain  inaccurately  located  peaks  or  suffer  some  peak  ] 

sfilitting  at  these  two  signal  loi'ations.  Whatever  occurs,  patU'rn  averaging  has  again  servinl 
to  maintain  [lattern  stability  so  that  signal  peaks  are  accuraU'ly  predicted  even  though  some- 
what broadened. 

In  a further  examination  of  the  effects  of  signal  phase,  another  averaged  MKSA 
antenna  (lattern  is  shown  in  Fig.  18  for  the  same  signals  and  noise  and  [larameters  of  the  two 
preceding  examples.  The  exce[)tion  is  that  the  signal  incident  at  60°  has  a relative  phase  of 
n ( 180°).  I'here  is  evidence  of  severe  peak  splitting  at  the  wider  angle  signals.  It  aiipears  that 
a signal  phase  of  w is  most  troublesome  for  MKSA  and  is  responsible  for  severe  loss  in  ac- 
curacy and  peak  distortion. 

Phase  distortion  effects  are  examined  further  for  two  closely  spaced  signals  separated 
by  6°  with  a S/rj  of  20  dB  each,  and  detected  with  an  8-element  antenna  using  a 6-point 
prediction  filU*r.  The  resulting  averaged  MKSA  antenna  pattern  for  the  two  signals,  both 
having  zero  initial  phase,  is  shown  in  Fig.  19  where  the  two  signals  are  separated  by  a shal- 
low null  of  only  a few  dB.  The  averaged  MKSA  pattern  for  one  signal  (the  one  at  +3°) 
having  an  initial  phase  of  ir!2  is  shown  in  Fig.  20.  Clearly,  the  effect  of  the  7r/2  relative  (and 
initial)  phase  is  to  shift  the  two  signal  peaks  further  apart  and  to  enhance  resolution  at  the 
loss  of  accuracy.  Also,  the  signal  peaks  appear  somewhat  broadened. 

In  Fig.  21  an  averaged  MKSA  pattern  is  shown  for  the  same  two  signals  in  the  two 
previous  examples,  but  with  the  one  signal  incident  at  +3°  having  an  initial  (and  relative) 
phase  of  tt  ( 180°).  The  F)attern  distortion  is  similar  to  that  depicted  in  Fig.  20  for  a relative 
phase  of  7t/2,  but  is  not  nearly  so  severe.  Distortion  effects  Eire  also  minor  (similar  to  Fig.  21) 
for  a fihase  shift  of  )t/4. 

Summary  of  Observed  MESA  Characteristics 

The  illustrations  of  MKSA  antenna  patterns  show  that  these  patterns  have  substantially 
reduced  [leak  widths  and  consequently  have  better  resolution  than  conventional  antenna 
jiatterns.  Also,  signal  detection  with  MKSA  appears  to  be  improved  over  signal  detection 
with  conventional  fiatterns.  Such  imjirovement  is  indicattxl  by  the  substantial  rwluction  in 
side  (noise)  peak  levels  achieved  with  postprocessor  integration  (pattern  averaging). 

Multiple  signals  of  varying  power  levels  may  be  accurately  predicted  with  MESA, 
although  the  number  of  antenna  elements  must  be  greater  than  twice  the  number  of  signals 
to  be  detect«>d. 
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Fig.  17  — Five  signals  (phase  at  60°  is  7T/4  ) 
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Fig.  19  — Two  signals  at  —3°  and  3° 
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Fig.  20  — Two  signals  (relative  phase  of  n/2) 
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Fig  21  — Two  signals  (relative  phase  of  n) 


Examples  of  phase  distortion  are  noted  for  phases  tt  and  rr/2,  which  are  thought  to  be 
most  troublesome. 

One  type  of  averaging,  postprocessor  integration,  is  shown  to  be  most  effective  in 
providing  stability  in  the  computed  MESA  patterns  and  in  substantially  improving  signal 
detection.  It  is  likely  that  preprocessor  integration  is  equally  effective,  and  may  be  prefer- 
able to  postprocessor  integration  when  computational  time  is  considered. 


OPTIMAL  FILTER  SIZE 

The  size  N (number  of  filter  coefficients)  of  a MESA  filter  is  constrained  to  be  at  least 
one  less  than  the  number  of  data  samples  M,  i.e., 

Af  < Af  - 1, 

so  that  at  least  one  data  sample  which  is  not  convolved  with  the  filter  coefficients  is  available 
for  estimating  the  prediction  error.  The  lower  Ijound  on  N is  dependent  upon  the  total  num- 
ber of  spectri  parameters,  since  some  minimal  number  of  filter  coefficients  is  required  to 
accurately  represent  all  spectral  component  parameters  such  as  amplitude,  frequency,  and 
phase.  For  instance,  if  there  are  P spectral  components,  all  with  the  same  relative  phase, 
then  N is  constrained  as  follows; 


2P  < Af  < A/  - 1, 


(26) 


where  2P  represents  the  total  number  of  spectral  parameters. 


KING 


While  N is  bounded,  the  actual  filter  size  is  optional  within  the  bounds  of  Eq.  (26). 
Anderson  [18]  and  others  have  noted  that  the  criteria  for  selecting  the  filter  size  depend 
upon  the  intended  application  or  function  of  the  MESA  power  spectra.  For  in  using  MESA, 
it  is  observed  that  the  spectral  characteristics  of  MESA  are  a function  of  the  filter  size.  Both 
resolution  and  peak  height  are  improved  at  the  larger  filter  sizes,  whereas  stability  and  ac- 
curacy sometimes  are  better  at  the  smaller  filter  sizes. 

Some  criteria  for  determining  the  filter  size  are  required  if  MESA  is  to  be  used  in  a 
completely  automated  manner  to  determine  unknown  power  spectra.  One  criterion, 
which  is  representative  of  two  MESA  characteristics,  is  the  output  S/t?.  Peak  height  and 
resolution  are  two  interrelated  properties  that  are  optimal  with  maximization  of  the  output 
S/t].  King  [20]  has  noted  that  maximization  of  the  output  S/rj  at  a spectral  peak  is  a reli- 
able criterion.  Another  criterion  developed  by  Akaike  [21]  has  been  investigated  by  Ulrych 
and  Bishop  [8]  and  found  to  be  only  partially  satisfactory. 

King  [20]  noted  that  the  output  S/tj  at  a spectral  peak  cOq  is  given  by 

(S/v)^  (27) 


where  the  power  spectra  are  evaluated  at  the  spectral  peak  a>o . With  use  of  the  power  spectra 
expression  (Eq.  6), 


(28) 


where 


1 -H 


N 

E 

n=l 


(29) 


A maximum  (•SA?)^  requires  that  Ewfcu^)  be  a minimum.  Therefore,  an  optimal  filter  size 
is  the  filter  size  (AT^)  that  minimizes  rj^/(cjQ). 

If  the  maximum  S/77  is  a criterion  for  determining  the  filter  size,  then  the  MESA 
properties  of  accuracy  and  stability  are  not  optimized,  and  remain  as  inherent  MESA 
characteristics.  Both  accuracy  and  stability  (under  varying  noise  fields)  have  been  satisfactory 
with  use  of  a maximum  S/77  output.  However,  when  signal  relative  phase  is  nonzero,  spectral 
peaks  of  the  computed  power  spectra  are  often  instable.  Sometimes,  averaging  of  such  com- 
puted spectra  appears  to  restore  the  stability.  Examples  of  MESA  antenna  patterns  of  phase- 
shifted  signals  are  provided  in  the  previous  section. 

Of  course,  maximization  of  the  output  S/77  at  each  spectral  peak  requires  that  the  peak 
locations  Wq  be  known.  The  peak  locations  may  be  determined  by  solving  for  the  roots  of 
the  function  r^(c>j)  at  a stable  and  accurate  low-order  filter  size.  Knowledge  of  the  peak 
locations  is  also  most  helpful  in  accurately  representing  computed  power  spectra. 
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An  example  of  the  value  of  maximizing  the  output  is  given  in  Fig.  22(a-g), 

where  the  MKSA  wavenumber  spectra  are  computed  for  a signaPlocated  at  0“  (broadside 
to  the  antenna).  The  antenna  has  eight  elements  and  the  input  S/tj  is  0 dB  per  element.  The 
MKSA  power  spectra  are  shown  computed  for  all  possible  filter  sizes  2 < N < 7 in  Fig. 
22(a-f). 

A familiar  problem  with  MKSA,  line  splitting  occurs  at  the  larger  filter  sizes  5 < A/  < 7. 
At  the  lower  filter  sizes  2 < ^ < 4,  the  jwwer  spectra  having  the  best  resolution  and  greatest 
jH*ak  height  occur  at  the  filter  size  N = 4,  Fig.  22(c).  The  filter  size  = 4 is  also  the  size  de- 
terminetl  by  maximizing  the  output  , and  the  optimal  power  spwtra  for  N = 4 are 

shown  in  Fig.  22(g).  Of  the  six  possible  filter  sizes,  the  power  spectra  for  N = 4 are  clearly 
those  with  the  best  resolution  and  peak  height.  Maximization  of  the  output  (S/r?)  also 
serves  to  avoid  such  problems  as  peak  splitting  as  observed  with  the  example  given.” 


CONCLUSIONS 

Kxamples  of  MKSA  anUmna  patterns  prest'nted  in  this  report  are  an  improvement  over 
conventional  phased-array  antenna  patterns.  Several  direct  comparisons  illustrate  the  im- 
proves! resolution  and  side-peak  suppression  that  are  characteristic  of  MKSA.  However,  it 
should  be  emphasized  that  these  results  were  obtained  from  computer  simulation  modelling 
of  the  antenna  elements  and  the  sources.  The  simulations  did  not  include  degrading  effects 
such  as  antenna  errors,  signal  scintillations,  etc.,  so  that  it  is  not  appropriatt*  to  assume  that 
th€*so  remarkable  resolution  characteristics  of  MKSA  will  hold  under  real-life  conditions. 

The  best  means  of  applying  MKSA  to  actual  radar  signals  has  not  been  established. 
Apparently  some  form  of  integration  (or  averaging)  is  necessary  at  some  point  within  the 
MKSA  algorithm  in  order  to  gain  stability  in  the  computed  antenna  patterns.  Postintegration 
(or  pattern  averaging)  demonstrated  in  this  report  is  one  method  for  obtaining  pattern 
stability.  However,  stability  may  also  be  obtained  by  averaging  the  elements  of  the  covari- 
ance matrix,  the  prediction  errors,  or  the  pn*diction  filter  coefficients. 

While  MKSA  appears  to  be  a useful  method  for  obtaining  improved  antenna  patterns, 
it  remains  to  qualify  these  improvements  and  to  establish  the  appropriate  applications  of 
MKSA  to  radar  problems  through  further  analytical  and  experimental  investigation. 
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APPENDIX  A 

WIENER  PREDICTION  FILTER 


PREDICTION  FILTER 

Consider  a signal  x(f)  with  a waveform  known  only  over  the  time  interval  (0,T).  The 
waveform  may  be  predicted  (or  estimated)  for  points  outside  the  window  (0,T)  by  using  the 
prediction  filter  in  a convolution  with  the  known  signal  x{t)  as  follows  (suggested  by  Wiener*) 

x(t)  = f a(T)x(t  - T)dT,  (Al) 


where  x(t)  is  the  predicted  signal  and  a(T)  is  the  impulse  response  of  the  prediction  filter. 
If  the  signal  has  been  adequately  sampled  within  the  time  interval  (0,T)  such  that 


< 2f, 


for  a sampling  interval  At,  M data  samples,  and  signal  frequency  f,  then  the  discrete  convolu- 
tion may  be  employed  as  follows; 


Xu  = 


N 


n=l 


N < M, 


(A2) 


where  the  summation  is  taken  over  N filter  coefficients,  N being  less  than  the  number  of 
data  samples. 


PREDICTION  ERROR  FILTER 

An  error  may  be  defined  for  the  known  set  of  data  samples  by 

N 

= ^ %^k-n  (A3) 

= Vfe-n 

n»0 


*N.  Wiener,  Extrapolation,  Interpolation,  and  Smoothing  of  Stationary  Time  Series,  John  Wiley  and  Sons, 
Inc.,  New  York,  N.Y.  1949. 


KING 


for 


7q  = l,n  = 0, 

where  7^  is  the  «th  prediction  error  filter  coefficient. 

The  squared  error  is  expressed  as  follows: 

N 

4-Z 

n=0 

Since  7q  = 1,  there  are  N remaining  unknown  prediction  error  coefficients.  These  N coeffi- 
cients may  be  determined  by  minimizing  the  total  mean  squared  error  which  is  defined  by 

M 

" M + 1 2 
*=-M 

The  summation  is  taken  over  all  errors  possible  to  compute  in  a forward  prediction  within 
the  data  window  defined  by  N + 1 data  samples.  The  incorporation  of  Eq.  (A4)  into  Eq. 

(A5)  results  in  the  following  equation; 


1 

N + 1 

M 

E 

N N 

"^m^k-n^k-m 

k=-M 

m=0  n~0 

N 

N 

M 

E 

E 

1 

M + 1 

^k  n4-m 

m=0 

n=0 

k=-M 

N N 

in 

11 

1 "S 

2 / 

C E 

n'^n^m  , 

m“0  n=0 


where  r_  represents  the  data  set  autocorrelation  coefficients. 

m’fi 

The  prediction  error  coefficients  are  defined  with  minimization  of  the  total  prediction 
error  as  follows: 

97* 

for  fe  = 1,  2,  3, . . . N.  The  resulting  N equations  with  N unknowns  are 


N 


n“d 


(A7) 


♦ 

^fe-m 


(A4) 
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fc  = l-To*-*  + Ti^o  + Tz*-!  + . . . + Tyv'-N  l = 0. 

* = 2,7o^2  T'l''*  ■"  >2''0  + • • • + yN'’N-2  ' 

k = + yx'-N  i * ')'2'’i-2  + • • ■ + 7jvro  = 0. 

An  additional  equation  which  defines  results  for  fe  = 0 where  Eq.  ( A6)  is  evaluated  as 
follows: 

N 

^N^yo'^yn^n  ioth=Q,  (A9) 

n*0 

when  it  is  recalled  that  7q  = 1. 

Bunit  has  assumed  that  the  mean  squared  error  power  E[^,  as  defined  by  Eq.  (A5) 
for  the  time  domain,  is  equivalent  to  the  mean  squared  error  power  which  is  given  by 
Eq.  (8)  for  the  frequency  domain.  However,  Kingt  has  noted  that  En  and  Pf^  are  equivalent 
only  when  the  predicted  noise  power  equals  the  actual  noise  power  or  when  the  noise  power 
is  negligible.  It  follows  that  for  high  signal-to-noise  conditions,  it  may  be  assumed  that 


Pn  = (AlO) 

and  Pn  may  be  evaluated  as  follows: 


N 


Pn 

= ^ 7^r^  for  fe  = 0, 

(All) 

n*0 

yo'^o  ^ yt 

+ Tz'-z  + • • . ys’-N  = Pn- 

(A12) 

If  Eq.  (A12)  is  combined  with  the  set  of  equations  (A8),  the  resulting  set  of  equations  may 
be  expressed  in  matrix  formulation  as  follows; 


''o''l'‘2''3  • • 

r 

>0 

Vor^rj  . . 

•'•n-1 

yt 

♦ ♦ 

fjrjrorj  . . 

• ’‘N-2 

♦ 

'•'2 

a IS  a 

'‘n’'nVn-2 

* 

'■yv  3 • • ■ '■q 

(A13) 


t'John  P.  Burg,  “Maximum  Entropy  Spectral  Analysis,"  presented  at  the  37th  Meet.  Soc.  Explor.  Geophys., 
Oklahoma  City,  Oklahoma,  Oct.  1967. 

tWilliam  R.  King,  "Some  Effects  of  Noise  Upon  Maximum  Entropy  Spectral  Analysis,”  NRL  Memorandum 
Report  3646,  Nov.  1977. 
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APPENDIX  B 
THE  BURG  TECHNIQUE 

Burnt  has  propost-d  a method  for  computing  a st't  of  N + 1 prtnliction  error  filter  coef- 
ficients as  a function  of  a known  si*t  of  N coefficients  using  the  Levinson  recursion  t>quation 


where  the  only  unknown  in  Kq.  (Bl)  is  the  last  coefficient  7^  + } of  the  new  set.  Burg  sug- 
gested that  the  unknown  coefficient  7^*  j be  obtained  in  a least-square  error  analysis  that 
incorporates  both  the  forward  and  backward  prediction  of  the  feth  point  as  follows: 


n=l 

N 


where  the  forwanl  prediction  and  the  backward  prediction  are  expresGi*d  as  discrete 
convolutions  of  the  prt*diction  filter  set  (of  size  N)  with  the  st't  of  A'  1 data  samples. 

It  follows  that  corresponding  forward  and  backward  prediction  errors,  denoted  by 
and  (jjJ'  respectively,  have  the  following  representations: 

«^  = X.  - X.  , 


N 


where  M is  the  number  of  data  samples,  and 

= X,  - X,  , 


N 


where  = 1 and  7^  = — . 


tjohn  P.  Burg,  "A  New  Analysis  Technique  for  Time  Series  Data,”  presented  at  NATO  Adv.  Study  Inst. 
Signal  Process,  with  Emphasis  on  Underwater  Acoust.,  Enschede,  Netherlands,  1968. 


34 


NRL  REPORT  8298 


forward  and  backward  prediction  errors  are  interrelated 
tnrougn  the  Levinson  recursion  relations 


= . 

+ 

N 

k ’ 

(Af  + 1 <Af), 

(B6) 

(^•0' 

. (/i 

+ (Jf  , 

, (l<;<Af-Af). 

{B7) 

computed  with  the  Levinson  recursion  relations  given  by  Eqs 
(B6)  and  B7)  m a bootstrap  manner  that  greatly  reduces  the  number  of  computations  re 
quired  with  use  of  the  matrix  formulation  given  by  Eq.  (10). 

The  total  mean  squared  error  may  be  expressed  as  the  sum  of  the  forward  and  back- 
ward mean  squared  errors  as  follows: 


E kJ  * k 


for  a filter  of  size  N where  N < M - 1,  with  Af  = number  data  samples. 

i!?  to  Eqs-  (B6)  and  (B7)  as  written,  the  total  mean  squared  error  may  be  ex- 
pressed for  filter  size  Af  + 1 as  follows.  ^ 


2(M-N~ 


M-N-if,  , .-I 

1 Y"  Y + 

N~l)  ^ \Pk  } 

fc*i  L>  / \ y_ 


i^g  e?prei£n  for  recursion  relations,  Eqs.  (B6)  and  (B7),  results  in  the  follow- 

/ \2  / ♦ \~ 

^iv-l  " 2(M--Ar-1)  ^ + ajn  (BIO) 


^ ALET(Y3VTRT5‘'6i''V^^^  f’*"*'  Maximum  Entropy  Spectrum  and  the  Burg  Technique  ’ 

ALEX  (03)-TR-75-01.  Texas  Instruments.  Inc..  Equipment  Gmup.  P.O.  Box  6015.  Dallas.  Texas  June  197^’ 
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In  order  to  minimize  the  mean  squared  error,  the  partial  derivative  with  respect  to  7^^  + } 
(the  only  unknown  filter  coefficient)  is  set  equal  to  zero: 


dt: 


N*\ 


= 0 


M-S-l 


2{M-N-  1) 


fc-l 


'^  + 1 nN 

Jn*\ 


k+N 


= 0.  (BID 


Kquation  (BID  niay  be  solved  for  7]^**  as  follows: 


(B12) 


The  new  larger  set  of  prediction  error  coefficients  may  be  evaluated  with  use  of  the  recursion 
relation,  Eq.  (Bl),  and  Eq.  (B12),  which  greatly  reduces  the  number  of  computations  re- 
quired of  the  matrix  formulation  given  by  Eq.  (10). 


